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Abstract. We study the low energy effective action for some collective modes of the 
color flavor locked phase of QCD. This phase of matter has long been known to be a 
superfluid because by picking a phase its order parameter breaks the quark-number 
U(1)b symmetry spontaneously. We consider the modes describing fluctuations in 
the magnitude of the condensate, namely the Higgs mode, and in the phase of the 
condensate, namely the Nambu-Goldstone (or Anderson-Bogoliubov) mode associated 
with the breaking of U(1)b- By employing as microscopic theory the Nambu-Jona 
Lasinio model, we reproduce known results for the Lagrangian of the Nambu-Goldstone 
field to the leading order in the chemical potential and extend such results evaluating 
corrections due to the gap parameter. Moreover, we determine the interaction terms 
between the Higgs and the Nambu-Goldstone field. This study paves the way for a 
more reliable study of various dissipative processes in rotating compact stars with a 
quark matter core in the color flavor locked phase. 
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1. Introduction 

At extremely high densities, Quantum chromodynamics (QCD) predicts that the 
individual nucleons that form standard hadronic matter should melt and their quark 
matter content should be liberated pQ. At the low temperatures expected in 
sufficiently old compact stars, quark matter is likely to be in one of the possible color 
superconducting phases, whose critical temperatures are generically of order tens of 
MeV [2], [3] . Since compact star temperatures are well below these critical temperatures, 
for many purposes the quark matter that may be found within compact stars can 
be approximated as having zero temperature, as we shall assume throughout. At 
asymptotic densities, where the up, down and strange quarks can be treated on an 
equal footing and effects due to the strange quark mass can be neglected, quark matter 
is in the color flavor locked (CFL) phase [HE]- The CFL condensate is antisymmetric 
in color and flavor indices and involves pairing between up, down and strange quarks. 
The order parameter breaks the quark-number U(1)b symmetry spontaneously, and the 
corresponding Nambu-Goldstone (NG) boson (Anderson-Bogoliubov mode) determines 
the superfluid properties of CFL quark matter. The gapless excitation corresponds to 
a phase oscillation about the mean field value of the gap parameter. The fluctuation in 
magnitude of the condensate is associated with a massive mode, which we shall refer to 
as the Higgs mode. Both of these fluctuations are excitations of several fermions and 
therefore describe collective modes of the system. 

The low energy properties of the system are completely determined by these 
collective modes [5], and their study is mandatory to understand dynamical properties 
that take place in compact stars with a CFL core. The actual superfluid property of 
the system is due to the presence of the gapless excitations, which satisfy the Landau's 
criterion for superfluidity [6] . Moreover, in rotating superfluids the interactions between 
NG bosons and vortices lead to the appearance of the so-called mutual friction force 
between the normal and the superfluid components of the system. The NG bosons 
are also responsible of many other properties of cold superfluid matter, in particular 
they contribute to the thermal conductivity and to the shear and bulk viscosities [6]. 
A study of the shear viscosity and of the standard bulk viscosity coefficient due to 
NG bosons has been done in Refs. [TJ [8] . A more detailed study of the bulk viscosity 
coefficients of the CFL phase has been done in [9] . In the CFL phase also the global 
SU(3)r x SU(3)l symmetry is spontaneously broken and the corresponding pseudo-NG 
bosons can contribute to the transport properties of the CFL phase. The contribution of 
kaons to the bulk and shear viscosity has been studied in Refs. |10l [TT] . The contribution 
of NG bosons to the thermal conductivity and cooling of compact stars were studied 
in Refs. [121 Ell EH US]- Other low energy degrees of freedom are the plasmons, which 
have been studied in Refs. [161 02] • 

Although the CFL phase is characterized by many low energy degrees of freedom, 
in the present paper we shall focus on the Nambu-Goldstone field, <fr, and on the Higgs 
field, p, because we aim to build the effective Lagrangian for the description of NG 
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boson-vortex interaction. Vortices are described by spatial variation of the condensate 
and therefore are determined by the space variation of the Higgs field. In principle one 
could study interaction of vortices with other low energy degrees of freedom, but the NG 
boson associated to the braking of U{1)b is the only massless mode, whereas pseudo-NG 
bosons associated to the breaking of chiral symmetry have a mass of the order of few 
keV [18], and plasmons associated to the breaking of SU(3) C symmetry have even larger 
masses [16], [T7]. Therefore all of these modes are thermally suppressed in cold compact 
stars. 

We derive the effective Lagrangian of the Nambu-Goldstone field, and of the Higgs 
field using as microscopic theory the Nambu-Jona Lasinio model [T9l 120] with a local 
four-Fermi interaction with the quantum numbers of one gluon exchange. This model 
mimics some aspects of QCD at large densities [2H [22] . We derive the interaction 
terms of the NG bosons up to terms of the type (<90) 4 . The leading contributions to 
these interaction terms were obtained by Son in Ref. [23] using symmetry arguments 
and the expression of the pressure of the CFL phase. To our knowledge the results 
of Ref. [23] so far have not been obtained starting from a microscopic theory. Only 
the free Lagrangian was obtained in [18]. We extend the results of [23] including the 
next to leading corrections proportional to the gap parameter. In order to do this we 
integrate out the fermionic degrees of freedom employing the High Density Effective 
Theory (HDET) [3]. Moreover we determine the kinetic Lagrangian for the Higgs field 
and the interaction terms between the Higgs and the NG bosons. 

A similar analysis to the one presented here was done for non-relativistic systems 
at unitarity in [23] and in [25] [26] . As we shall discuss in Section [4] the main difference 
with the non-relativistic systems at unitarity is that in the CFL phase integrating out 
the Higgs field does not lead to a change of the speed of sound. 

This paper is organized as follows. In Section [2] we present the Nambu-Jona Lasinio 
model and we determine the expression of the effective action of the system in the HDET 
approximation. In Section [3] we derive the effective Lagrangian for the NG boson, 
neglecting the oscillations in the magnitude of the condensate. In Section H] we derive 
the LO interaction terms between the Higgs mode and the NG bosons and the kinetic 
Lagrangian for the Higgs field. We also compare our results with the non-relativistic 
results of [23] . We draw our conclusions in Section [5] Some interaction vertices as well 
as the calculation of some integrals are reported in the Appendix. 

2. The model 

We consider a Nambu-Jona Lasinio (NJL) model of quark matter with a local four- Fermi 
interaction as a model of QCD at large quark chemical potential fi. We assume that 
fi ^> m s and therefore we neglect u, d and s quark masses. In the absence of interactions 
the Lagrangian density describing the system of up, down and strange quarks is given 
by 
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where i,j = 1, 2, 3 are flavor indices, a, (3 = 1, 2, 3 are color indices and the Dirac indices 
have been suppressed. In QCD one can show that at large densities the gluon exchange 
leads to the formation of a quark-quark condensate [2]. At asymptotic densities the 
favored phase is the color flavor locked phase [I] which is characterized by the condensate 

<v&c^>~a £ e ^' 7 , (2) 

1=1,2,3 

which locks together color and flavor rotations. Here C = ry 2 7° is the charge conjugation 
matrix and e a pi and e 1 ^ 1 are Levi-Civita tensors. This condensate breaks several 
symmetries of QCD, see p[| for a detailed analysis. For our purposes it is sufficient 
to note that the CFL condensate breaks the U(1)b symmetry and therefore leads to the 
appearance of a gapless Nambu-Goldstone boson. 

In the NJL model the interaction among quarks mediated by gluons is replaced by 
a four-Fermi interaction of the BCS type. We shall consider an interaction with the 
same quantum numbers of one gluon exchange 

Ci = ~g^ A Hh^X A ^, (3) 

where g > is the coupling constant and \ A with A = 1,...,8 are the Gell-Mann 
matrices. This interaction at large chemical potentials leads to the formation of the 
CFL condensate. Introducing the new basis for the quark fields 

1 9 
V2 A=1 



where A = y 2/3 x I, we have that 

{^ A C^ B ) ~ A AB , (5) 

where A A b = A A 5 AB , where Ai = ... = Ag = A and Ag = — 2A. In this basis the 
four-Fermi interaction is given by 

A = —VABCDeab^a^b ^f^T > ( 6 ) 
8 

where Vabcd — Tr '^2(\a^e^b^c^e^d) an d a(a) = 1,2 are the Weyl indices for L(R) 

E=l 

components [2Tj . 

In order to study the fluctuation of the condensate we introduce the Hubbard- 
Stratonovich fields A A b{x) and A* AB (x) which allow to write the partition function 
(normalized at the free case for g = 0) as 

1 ■ J [dA, dA*} exp 1 1 j d 4 x 



A A bW AB CD^cd , r 
V L-A 
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where 

W A bcdVcdef = o~aeo~bf and VabcdWcdef = o~aeo~bf , (8) 
and the semi-bosonized Lagrangian is given by 

C A = M$ + /^7o)^a - \^ABi^ A Cr B ) + I&ab&aWb) ■ (9) 



(7) 
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The fluctuations of the condensate A(x) around the mean field value A MF can be 
described by two real fields; one field describes the variation of |A(x)| while the other 
field is associated to a local phase change. Therefore we write 

A AB (x) = IA%£ + p AB (x)]e 2 ^ , (10) 

where pab(x) = p(x) diag(l, 1, —2)ab and 0(x) are the real fields. Hereafter we shall 
suppress the indices A, B and we shall indicate with A the mean field value of the gap 
in order to simplify the notation. We find convenient to redefine the fermionic fields as 

^^^e^ (x) , (11) 
and in this way the semi-bosonized Lagrangian is given by 
C A = $ (V<9 M + tV " 7°<9 O - 7^0) i> 

- ~^C(A + p)^ + \^C(A + p)j, . (12) 

Writing the Lagrangian in this form it is possible to define an "effective chemical 
potential" 

fi = p-d (p. (13) 

In other words do4> describes long-wavelength fluctuations of the chemical potential 
on the top of its constant value /i. We want now to clarify one point regarding the 
terminology used for the field <fi. Since d$(f) describes fluctuations of the chemical 
potential, it is not correct to call this field the phonon, which describes pressure 
oscillations. However, when one neglects the effect of the gap A, the Lagrangian of 
the phonon and of the NG boson associated with the breaking of U(1)b symmetry 
coincide. The reason is that for vanishing values of A the oscillations of pressure are 
proportional to the oscillations of the chemical potential. 

It is convenient to define a fictitious gauge field = (<9o0, V0) which according 
with Eq. (TT2|) is minimally coupled to the quark fields. Hereafter we shall refer to A^ as 
the gauge field, although it is not related to any gauge symmetry of the system. With 
this substitution we rewrite the Lagrangian in Eq. (TL2l as 



C A = $ (i-fDp + /i 7 °) $ - \^C{A + p)r + \^C{A + p),f, , (14) 



where 0^ = 8^ + iA^. 



In order to simplify the calculation we employ the High Density Effective Theory 
(HDET), see [3]. Using standard techniques, the Lagrangian describing the kinetic 
terms and the interaction with the gauge field A 11 can be written as 



dw 

8^ 



t / P^ U D U D 
iplliV-D ^ 




+ rjA (iV-D- P ^ D » D » ) ^_ (15) 
^~ V 2p + iV-D> 

where the positive energy fields with "positive" and "negative" velocities are given by 
^ ± = M±v), (16) 
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and where 

^ = (l,v), V"* = (l,-v), P^ = g^ 1 , (17) 

with v the Fermi velocity. The non-local interactions which are present in Eq. (1151) are 
due to the integration of the negative energy fields [3]. 

In order to have a compact notation it is useful to introduce the Nambu-Gorkov 
spinor 

* = — ( ^ + \ (18) 

which allows to write the semi-bosonized Lagrangian as 

U = C a + C b , (19) 

where 

rdv ,(iV^-V»A» A + p \ 

is the sum of the kinetic term and of the local interaction term, while the non-local 
interaction term is given by 

where L = (2/i + iV ■ D)(—2/i + iV ■ D*) — A 2 — ie. More details about the derivation 
of the non-local interaction will be given in [28] . 

The semi-bosonized Lagrangian is quadratic in the fermionic fields and therefore we 
can integrate them out. In this way the effective action can be written in terms of the 
fields <p and p and their derivatives. An analogous calculation in the non-relativistic case 
has been done in [25] . Before doing this we consider in more detail the interaction terms 
between the fermionic fields and the fictitious gauge field determined by the non-local 
term. In momentum space we have that 

C b = f d ^P^Aj- 2 " + V -' +V - A ~" A - U*,(22) 

b J 4vr M "\ -A 2fi + V-£-V-AJL ,KJ 

where in momentum space L = (2/i + V ■ i — V ■ A)(—2fi + V ■ £ + V ■ A) — A 2 — ie. The 
"residual momentum" of the quarks, i^, is defined as follows: 

4 = Po £i=Pi~ IMi , (23) 

where p^ is the four-momentum of the quarks. 

Expanding the denominator of the expression in Eq. ( 1221) in powers of A and 
considering terms up to the order A 4 , we have that 

Cb = I £ {^ V 2 V ^A,A V + &r$">*A,A v A p 

+ ^T^ pa ^A^A u A p A a )+0(A 5 ), (24) 
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Figure 1. The tree- level interaction vertices of the gauge field with the fcrmions in the 
HDET. Full lines correspond to fermionic fields; wavy lines correspond to the gauge 
fields. The vertices are named Ti, T2, T3 and Fi, where the subscripts indicates the 
number of the external gauge fields. Their expression is reported in the | Appendix A| 



Figure 2. Interaction of the Higgs field (dashed line) with quarks (full line). The 
corresponding vertex is named T p and is reported in Eq. (|29[) . 



where the expression of the vertices T 2 , T 3 and T 4 are reported in the Appendix A We 
also define 

r -(T A)- <*> 

which is the vertex that describes the minimal coupling of quarks with the gauge field 
in the HDET. The various vertices are schematically depicted in Fig. [TJ with the wavy 
lines corresponding to the gauge field and the full line corresponding to the fermionic 
field. 

Integrating out the fermionic fields, the partition function turns out to be given by 
Z /[<<A,dA-]exp[i/d 4 x A.ibW^bctAU det[S-'] 1/2 

% = dSi^F s exp [i51 (26) 

where S is the action of the system and the full inverse propagator is given by 

S- 1 = S- M \ + r . (27) 
The mean field inverse propagator is given by 

S ™ F = ( 1 V -A iV%, ) ' (28) 

while T = T p + Yi + T 2 + T 3 + r 4 . The interaction of quarks with the p field is described 
by the vertex 

r -=(-,7)- < 29 ' 

while the interaction of quarks with the NG bosons is given by Ti = T^A^, T 2 = 
ApAj,, T 3 = r 3 p A^A u Ap and T4 = T^ p A^A v ApA a . 
We separate the mean field action from the fluctuation writing 

S = Smf + <S e s , (30) 
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where <S e g is the effective action describing the low energy properties of the system. The 
mean field action provides the free energy of the system 

n = S MF = - l -Tr HSmf] - ~ [A AB W ABCD A* CD } , (31) 
^ 9 

and is a function of the quark gap parameter A. Hereafter, Tr, symbolizes the trace 
over the Nambu-Gorkov index, the trace over color-flavor indices, the trace over spinorial 
indices and the trace over a complete set of functions in space-time. 

The gap parameter can be determined by the stationary condition of the mean field 
action 

dA A~ ' 1 ] 

which in turn allows to determine the pressure of the system by P = —Q(A). In the NJL 
model, the mean field value of the gap parameter turns out to depend on the coupling 
g and on the three-momentum cutoff, see e.g. [5]. Since we do not need the numerical 
value of the gap, we shall treat A as a free parameter. 

The fluctuation around the mean field solution are described by the effective action 

Scff = - - / d 4 x [pab(x)W ABC dPcd(x) + 2p ab (x)W AB cdA C d} 
9 J 

- ^Trln(l + £Wr) , (33) 



which contains the T expansion 

Trln(l + 5 A/F r) = Tr 



00 (—~\) n + 1 



(34) 



and we shall evaluate terms up to the fourth order in the gauge fields. We shall also 
determine the leading order terms of the p field Lagrangian and the leading order 
interaction terms of the p field with the NG bosons. In the expansion we shall neglect 
the NLO terms of the kind (dd<p) n where n is any nonzero integer. 

3. The Lagrangian for the NG boson 

Neglecting oscillations in the modulus of the condensate we can determine from the 
Eqs. ( 133]) and ( |34|) the effective Lagrangian for the NG bosons. As we shall show below, 
for vanishing values of A we obtain the same results obtained in Ref. [23] 

^ = ^[(/i-9 o 0) 2 -(^) 2 ] 2 . (35) 

This expression relies on symmetry considerations, in particular on conformal symmetry, 
and on the expression of the pressure in the CFL phase. We shall reproduce these results 
and evaluate the corrections of the order (A//i) 2 which are related with the breaking of 
conformal symmetry. Our strategy is to first expand the Lagrangian in the fields, 
and write 

= A + £ 2 + £ 3 + A , (36) 
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where C m is the term with m gauge fields. Then we expand the various terms in A//x. 
Notice that the term C m will be obtained by the expansion in Eq. (134)) considering all 
the terms with n < m. 

A different way of obtaining the leading order in \x Lagrangian for terms 
proportional to do<fr is the following. At the leading order in \i — neglecting terms 
proportional to A - - the free-energy density of the CFL phase in the absence of 
oscillations is given by 

3 A 



n 



Air 2 



(37) 



We have seen in Section |2]that <9o0 corresponds to a fluctuation of the chemical potential 
of the system. Therefore, including these fluctuations the free-energy density of the CFL 
phase at the leading order in /i is given by 

(3* 



n 



4vr 2 



H 



where /i is defined in Eq. (fTBj) . Expanding the free energy we obtain that the Lagrangian 
for the 8q4> field is given by 



C(d 



47T 2 



7T 



9 /i 2 W) 2 -^(9o0) 3 



2tt 2 



7T 



471-2 



(a o 0) 4 -(39) 



3. 1 . One-point function 

The evaluation of the term of the effective Lagrangian proportional to dip requires the 
computation of the diagram in Fig. |3j One has that 

& = -iTr[ST]\ A , (40) 

where the subscript means that in the evaluation of the trace only the terms linear in 
A must be included. We find that 



where we have defined 



V ■ A- 



V ■ 



t 



dC 



7T 



E 

N c ,N f 



D* 

cW r+ 5 dh 



V ■ A- 



V ■ 



D 



An J~8 2ix 



°° Mo 

-oo 2tt 



(41) 



(42) 



where 5 is a cutoff that we shall set equal to /i, see [3]. In Eq. (14 ip we have also used [29] 
D = V ■ £V ■ £ - A 2 + ie , (43) 
and for future convenience we define the quantity 



L = (2/i + V ■ £)(-2fi + V ■ I) - A 2 



ie . 



(44) 



The sum over flavor and color degrees of freedom is straightforward and we obtain 



2i 

IT 

- V-A 



dv 

Air 

„ V ■ 



d 2 



(27T) 



V ■ A 8 



V ■ 



+ 



V ■ 



D(A)* D(-2Ay 



V ■ 



D(A) D(-2A) j 



(45) 
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Figure 3. One loop diagram that contributes to the coefficient of the do<j> term in the 
effective Lagrangian C\ in Eq. (|36|) . The full line corresponds to the quark field. The 
wavy line corresponds to the external gauge field. 



Employing the expressions reported in the Appendix B one can do the integration over 
the residual momentum and on the Fermi velocity. In this way one finds that 



3 n 6 9 . 

-— /i 3 + — /i 2 A 



<9 n 



(46) 



7T^ 7T Z 

Notice that at the leading order in fi the expression above for the term proportional to 
do<f) agrees with the corresponding expression reported in Eq. (|39|) . 



3.2. Two -point function 

The Lagrangian involving two gauge fields is given by 

A 2 



C 2 = -i[ Tr[ST] - -Tr[SrSr] 



(47) 



where in the evaluation of the trace one has to consider only terms quadratic in A. The 
corresponding diagrams are reported in Fig. HI The diagram in Fig. Hh) gives 

'A 2 + V-i(V -£-2/2) 



^[S?]\ A2 = ^dlP^A,A v 



LnD 



+ (V -)■ V) 



(4? 



where Lq has been defined in Eq. ( )44|) and (V — > V) actually means an expression that 
is obtained by replacing (V — > V, £\\ —> —£\\, e — > — e). Hereafter we shall always use this 
way of writing in order to simplify the notation. 





a) 



b) 



Figure 4. One loop diagrams that contribute to the coefficients of the Aq term, and 
of the A 2 term of the effective Lagrangian £2 in Eq. (|36|) . Full lines are fermionic 
fields. Wavy lines are gauge fields. 



The diagram in Fig. Hb) gives 



Tr[STST]\ A2 =f: 



\V-A)\l + £ ll ) 2 (V- AV ■ A) A' 



U 1 



+ (V -> V) 



(49) 
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and evaluating the integrals using the expressions reported in the Appendix B we have 
that 



2tt 2 
1 



1 - 2- 



A 5 



A 



2tt 2 



1-2.1^ A 2 + o(^log(A//i) 







:m D A 



m 2 M A 2 . 



Then the speed of the NG boson is given by c s = vtlmI^d and the corresponding plot 
is reported in Fig. [5j For vanishing A we have that c s — 1/ a/3 meaning that the system 
is scale invariant. The effect of a nonvanishing A is to increase the speed of the NG 
boson. 



0,36 



0,34 



0,02 0,04 0,06 0,08 0,1 
A/JLL 



Figure 5. Speed of the NG boson associated to the breaking of U(1)b as a function 

Of A/fJL. 



3.3. Three-point function 



The Lagrangian describing the interaction of three NG bosons is formally given by 



C, 



Ti[ST] - -Tr [5T5T] + iufSrSTST] 



A3 



(51) 



and corresponds to the evaluation of the diagrams in Fig. [6J The diagrams Ek) and 
[6b) contribute exclusively to the term A A 2 , while the diagram [6b) gives the term 
proportional to Aq. The contribution of the diagram in Fig. [6^) is given by 

~L V -£+(V-£ — V-£ — Ay) (D + 2A 2 - 2yV 



Tr l ST ]\ A3 = ^diP^A^V-A 



DL 2 



- (V^V) . 
The contribution of the diagram in Fig. [6b) is given by 

"2A 2 t> -£+(V -£-2y)((V 



(52) 



Trfsrsr] 



A 3 



2 } di P^AuAyV ■ A 



-A 2 



D 2 Lc 
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a) 



b) 



c) 



Figure 6. One loop diagrams that contribute to the effective Lagrangian £3 in 
Eq. ([55)) . Notice that the only contribution to Aq comes from diagram c). Full lines 
are fermionic fields. Wavy lines are gauge fields. 



- (V V) . 
The contribution of the diagram in Fig. [6b) is given by 

-{V ■ Af{V ■ if - A 2 (V ■ AfV ■ AV 



Ti[STSTSY] 



A 3 



P 



D ? 



+ 



2A 2 (V ■ A) 2 (V ■ A)(V 
LP 



(54) 



Evaluating the integrals with the help of the expressions reported in the Appendix B 
one finds that 

3A 2 ^ 



C 



3 ~~ ~~ 2 

7T 



2fi 



A 6 
2 l A) ■ 



(55) 



3.4- Four-point function 



The Lagrangian describing the interaction of four NG bosons is formally given by 
£ 4 = -i (Tr[ST] 



-Tr[srsr] + -Trf^rsrsr] + -i^srsrsrsr] 



2 3 1 4 ■ 'J A * 

and corresponds to the sum of the diagrams reported in Fig. [3 All these diagrams, with 
the exception of the diagram [7k-), originate from the non-local vertices r 2 , T 3 and T 4 , 
and therefore can give contributions to the coefficients of the terms with AqA 2 and A 4 . 
The diagram in Fig. [7b) contributes to the coefficient of the term proportional to A 4 . 
The contribution of the diagram in Fig. [TJi) is given by 



Tv[ST] 



A* 



+ z 



MP^A^ V (V-A)UV-A) 
(V • I — 2fi) 2 + (V ■ 1 + 2u) 2 



-V -£(V -l + 2u) 



DLl 



+ V ■ A 



2Z-A 2 



2Z 



(V ■ £ - 2fi){V ■ £ + 2y) 



+ (V -> V)\ ,(57) 



where Z = (D + 2 A 2 - 2fiV ■ 
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a) 



b) 



c) 





d) 



e) 



Figure 7. One loop diagrams that contribute to the effective Lagrangian £4 in 
Eq. (|36|) . Full lines are fermionic fields. Wavy lines are gauge fields. 



The diagram in FigJTb) gives 



dlP^A.AAiV -Af 



2A 



2 (V ■ if + 3(V ■ if - ZpV ■ i + 5fjV -£ + 4fj? 



D 2 L\ 



+ 2V ■ AV ■ A 

+ (V^ V)} , 
while the diagram in FigJTfc) gives 



diP^ApAvPoPAaAp 



(D + 2A 2 ) 2 + n(y ■ i(V ■ if - A 2 {V -i-3V-i))- 4(V ■ if/i 2 ' 



D 2 L 2 



(58) 



D + 2A 2 + 4// 2 (y 



4A 2 /i(2^ • i + n) 



D 2 L 2 



+ 



V ■ i(6A 2 V ■ £ - 4fi(V • if + AA 2 ^) 



D 2 L 2 



+ (v -»• y) 



Both these diagrams are determined by evaluating 



Tr[srsr] 



.A 4 



(59) 
(60) 



The diagram in FigJTU) gives 
TrfSTSTSr] 



+ A- 



MP^A^Uy ■ Af 
2{V ■ if + 7(V -if - 



3{V ■ if(V ■ i - 2a) 
D*Lo 
2fiV ■ £ + AfiV ■ i 



D 3 L 

9D + 12A 2 + Afi(V ■ i - 2V ■ 



V ■ AV ■ A 



DLl 



(61) 



Finally, the diagram in FigJT^) gives 

TrfSTSrSrST]^ = idijyi [{V ■ Af(V ■ if - 3(V ■ A) 3 V ■ AA 2 {V ■ if 

- V ■ A(V ■ AfA 2 (V ■ if + (V ■ Af(V ■ Af(2D + 3A 2 ) 
+ (V -> V)] . (62) 
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After a long but straightforward calculation we find that at the leading order in p 

In this case we restricted our calculation, for simplicity, to evaluate the leading terms 
in p. 



4. Higgs field and interaction terms 

The fluctuations in the absolute value of the condensate are described by the Higgs field 
p defined in Eq. f lTUj) . The effective Lagrangian for this field, including the interaction 
terms with the NG bosons, can be determined by means of the same strategy employed 
in the previous section. 




Figure 8. Diagrams that contribute to the one loop self-energy of the Higgs field. 
Full lines are fermionic fields. Dashed lines correspond to the Higgs field, p. 



From the V expansion in Eq. ( 1331) we obtain the self-energy diagram in Fig. [8] which 
gives for vanishing external momentum of the Higgs field 

C pp = -^p 2 p{x)\ (64) 

which is the mass term for the p field. However, the actual mass of the Higgs is not 
proportional to the chemical potential: as we shall show below one needs a wave function 
renormalization in order to put the Lagrangian in the canonical form. We notice that 
the p field does not carry color and flavor indices. In order to obtain Eq. (I64p (and all 
the expressions below), color and flavor indices have been properly contracted and the 
resulting Lagrangian is expressed in terms of colorless fields. 

Considering the fluctuation of the chemical potential given by do<f), the Lagrangian 
above turns into 



C pp = -- 2 ~p 2 p{xf = — (ai - <9 O 0) 2 P(*) 2 , 



which automatically gives part of the interaction of the Higgs field with the NG bosons. 
The remaining interaction terms of two Higgs field with two NG bosons correspond 
to a coupling (di<fr) 2 p(x) 2 . These interaction terms, as well as the interaction terms in 
Eq. (]65p . can be obtained from the diagrams in Fig. |HJ 
Evaluating these diagrams one obtains that 



> ~'Pp4>4> o 

7T Z 



{p-d^f-^f 



p{xf , (66) 



which describes the interaction terms among the Higgs fields and two NG bosons at the 
leading order in p. Notice that the term pdu(j)d^(j) describing the interaction between 
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Figure 9. Interaction of the Higgs field with two NG bosons at the one loop level. 
Full lines correspond to fermionic fields, dashed lines correspond to the Higgs field and 
wavy lines correspond to the gauge field = <9 M 0. 



one Higgs field and two NG bosons is missing. This is basically due to the fact that 
one cannot have a term like p?p from the loop expansion. The fact that in our case 
terms linear in the p field are missing, leads to an interesting effect. Let us consider a 
Lagrangian with terms up to quadratic order in the NG bosons and Higgs fields. Since 
the Lagrangian is quadratic in the massive Higgs field one can integrate them out from 
the theory obtaining the actual low energy Lagrangian of the system. In the CFL phase 
this does not lead to a modification of the effective Lagrangian of the NG bosons. In 
particular the velocity of the NG bosons at the leading order in \x is c s = This is 
rather different from what happens in the non-relativist ic case [21], where it was shown 
that integrating out the Higgs mode one obtains the modification of the speed of sound 
first determined in [30] . 

In order to clarify this point, we compare our Lagrangian with the quadratic non- 
relativistic Lagrangian. Schematically the result of [21] can be written as 



Cn.r. (p,Y) = |p 2 + BpY + ^Y 2 + DY , (67) 

where A,B,C,D are some coefficients and Y = d^cj) + , with m the mass of the 
non-relativistic fermions. The speed of the NG bosons is given by 

Integrating out the p field, the effective Lagrangian for the field turns out to be 

B 2 + AC n 

Cn.r. (0) = 2A Y + DY , (69) 

and the speed of sound is modified to 
2 _ DA 

° s ~ ~m{B* + AC) ' ( } 

Notice that if the coupling B vanishes, then the speed of the NG bosons remains the 
same one has in Eq. fl68l) . 

In our case the LO Lagrangian quadratic in the p field is given by Eq. fl66l) and there 
is no term that couples one p field with two NG bosons, i. e. expanding the Lagrangian 
one finds that the analogous of the coefficient B in Eq. (|67|) is missing. Therefore, 
integrating out the p field does not change the effective Lagrangian for the NG bosons. 
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From the diagrams in Fig. H] we can determine the kinetic terms of the effective 
Lagrangian of the Higgs filed. Considering soft momenta of the Higgs field, p << A, 
and expanding up to the order (p/A) 2 , we obtain 



13p z 1 



(doP? -\{d iP f 



1 12// 



2 



7T 2 



P 2 • (71) 



2 4tt 2 A 2 

By a wave function renormalization we can cast the above expression into canonical 
form and we obtain that the mass of the Higgs field is given by 

m p = 4A, (72) 

which means that the mass of the Higgs is twice the fermionic excitation energy. This 
is analogous to the result obtained in the chiral sector, where the masses of the mesons 
turn to be equal to twice the effective mass of the quarks, see [20] . 

From the expression of the mass of the Higgs, it seems unlikely that an NG boson 
could excite the p mode. The reason is that the scale of the field is T and in compact 
stars T < A. Even if thermal NG bosons in general cannot excite the Higgs, for vortex- 
NG boson interaction one has to properly take into account the fact that as one moves 
inside a vortex, the actual value of the mass of the Higgs field should decrease. The 
reason is that as one moves inside a vortex, the value of the condensate becomes smaller 
and smaller and the mass of the p field should decrease accordingly. Therefore at a 
certain point it should happen that m p < T and it may become possible for an NG 
boson to excite the p mode. However, a vortex is a modulation of the p field itself, 
and therefore the discussion of the interaction between NG bosons and Higgs field in a 
vortex is subtle. We postpone to future work the analysis of this situation. 

We now try to a give a general expression of the interaction terms between two NG 
bosons and the Higgs fields. Let us first neglect space variations of the field. Then, 
expanding the effective action one has terms like 

£{p, d <f)) = ^ c n — — , (73) 

n>2 

where c n are some dimensionless coefficients. This expression is simply due to the fact 
that any term in the effective Lagrangian is multiplied by p 2 , which comes from the 
phase space integration, while the denominator comes from dimensional analysis. 

Now we want to derive the expression of the terms with space derivatives of the 
NG bosons. We know from Eq. (166j) . that for n = 2, one has to replace (p — do<fr) 2 with 
(p — d (p) 2 — l/3(V</>) 2 , where the coefficient 1/3 is precisely the square of the speed 
of sound. In other words the correct metric for the propagation of NG bosons is the 
acoustic metric, = diag(l, — 1/3, — 1/3, — 1/3). Then, we introduce the four vector 
X M = (p — do<f), V</>) and our guess is that the Lagrangian at the leading order in p is 
given by 

Y u n n n 

C( P ,xn = E-r X X lf ■ (74) 

n>2 

The reasoning is that any perturbation propagating in the medium feels the presence 
of the background which induces the acoustic metric g pu [3T], |32]- We shall investigate 
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in more detail this guess in future work, however it was already shown in [52] that the 
free Lagrangian of NG bosons in the CFL phase can be written employing the acoustic 
metric. Therefore Eq. ( 174|) seems to be an educated guess. 

If the expression above is correct it allows to readily determine the interaction 
between any number of Higgs fields with two NG bosons in a straightforward way. As 
an example we can determine the interaction between three Higgs field and two NG 
bosons. We first evaluate the diagram in Fig. (TO], which gives the interaction among 
three Higgs field 

then, upon replacing fi 2 — > X lx X u g llv we obtain the LO interaction Lagrangian 

r 2X ^X V (jav 3 / — 

ZpppW = P ■ ( 76 ) 

It is not clear whether this reasoning can be extended in order to include subleading 
corrections of the order A//x. In principle one would expect that the leading effect should 
be to perturb the metric g^, but further investigation in this direction is needed. 




Figure 10. Interaction of three Higgs fields. Full lines correspond to fermionic fields, 
dashed lines correspond to the Higgs fields. 



5. Conclusions 

The low energy properties of cold and dense color flavor locked quark matter are 
determined by the NG bosons associated with the breaking of the U(1)b symmetry. 
At the very low temperatures expected in compact stars, the NG bosons probably 
give the leading contribution to the transport properties of the system. Therefore the 
detailed knowledge of their effective Lagrangian is important to precisely determine the 
transport coefficients. We have determined the effective Lagrangian for this field starting 
from a microscopic theory. As high energy theory we have considered the Nambu-Jona 
Lasinio model with a local four-Fermi interaction with the quantum numbers of one 
gluon exchange. Then we have gauged the U{1)b symmetry introducing a fictitious 
gauge field. Finally we have integrated out the fermionic degrees of freedom by means 
of the HDET. 
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We have confirmed the results of Ref. [23J, which were based on symmetry 
arguments, and extended including next to leading terms of order (A//i) 2 . These 
corrections are relevant because are related with the breaking of the conformal symmetry 
of the system and can be significant large for matter at non-asymptotic densities. For 
this reason, the present work paves the way for a more detailed calculation of the 
transport properties of CFL quark matter, including the scale breaking effects. 

We have also determined the interaction of the NG bosons with the Higgs mode, 
i.e. with the collective mode associated with fluctuations of |A|. These interactions are 
relevant in the calculation of the interaction of the NG bosons with vortices. Indeed, 
vortices can be described as space modulation of | A|. A preliminary study of the mutual 
friction force in the CFL phase has been done in Ref. [33] . However, in that calculation 
of the mutual friction force only the elastic scattering of NG bosons on vortices has 
been taken into account. Since we have determined the full low energy Lagrangian, 
we are now in a position to evaluate in more detail the interaction of NG bosons with 
vortices, including non-elastic scattering, however we postpone the analysis of vortex- 
phonon interaction to future work. Indeed the treatment of this interaction is non-trivial, 
moreover superfluid vortices which wind the U{1)b are topologically stable but according 
to the result of the Ginzburg-Landau analysis of Ref. [M] they are dynamically unstable. 
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Appendix A. Interaction vertices 

In the HDET fermions interact with gauge field by the following vertices. The minimal 
coupling is due to the vertex 



Non minimal couplings are due to the expansion of the non-local interaction in Eq. ( 12~TT) . 
Expanding in the number of gauge fields one has the term with two gauge fields: 




(A.l) 




(A.2) 
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three gauge fields: 

r 3 = ~W ■ A(2/2 + V - l) + V ■ A(2fi -V-£)]- ^P^A^ , (A.3) 
and finally the interaction with four gauge fields: 

r 4 = ^-[V ■ A(2fx + V ■ £) + V ■ A(2n - V ■ £)] . (A.4) 
In the present article we do not consider interactions with more than four gauge fields. 

Appendix B. Integrals 

In this appendix we evaluate some integrals used in the evaluation of the effective 
Lagrangians for the NG boson and for the Higgs field. For the integrals in £ and £u we 
define 



We have that 



d 2 £ = I d£ tt I d£ . 



d 2 £ * -i = (-l) n+1 ^- forn > 1 
D(£) n+1 v ; nA 2n 

d z £ ,J = / d z r ,J = -in S n i , 
D{£) n J D{£) n 



where D(£) = V ■ £V ■ £ — A 2 + ie. We also need the integrals 



L (£) 2 2 A 2 \ a/^ 2 + A 2 v/9/i 2 + A 2 J ' 
2 1 _. fi ( 2/i 2 + 3A 2 9(6u 2 + A 2 ) 

CL K. , „. — Z7T - 



where 



L (£) 3 8A 2 \(fi 2 + A 2 ) 3 / 2 (9/i 2 + A 2 ) 3 / 2 



L = (2/i + V ■ £){-2n + V ■ £) - A 2 - ie . 



B.l) 

B.2) 
B.3) 

B.4) 
B.5) 
B.6) 

B.7) 



Appendix B.l. Angular integrals 



Considering a general vector A^, for the integrals involving terms of order A we have 
that 

H = J^-V-A = A . (B.8) 
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At the order A 2 we have 

dv „,.„ . 



Hi= J -^P^A,A U = --K\ (B.9) 
H 2 = J d ^-{V.Af = Al + \A\ (B.10) 

tf 3 = J d ^V-AV-A = Al- l -A\ (B.ll) 

At the order A 3 we have 

r dv 9 
H 4 = J — P^A fl A l/ V-A = --A A 2 , (B.12) 

H 5 = J d ^{V-Af = Al + A,A\ (B.13) 

H,= J d ^{V-A){V-Af = Al- l -A Q A 2 . (B.14) 
At the order A A we have 

H 7 = J ^(V ■ A) 4 =A* + 2A 2 A 2 + ^A 4 , (B.15) 

H 8 = J ^-(V -A)(V- A) 3 = At- \a 2 A 2 , (B.16) 

H 9 =J ^(V ■ A)\V ■ Af = Aq- 2 -A 2 A 2 + iA 4 , (B.17) 

r dv ~ 9 9 

H w = J ^P^A,A U (V ■ A) 2 = --A 2 A 2 - ^A 4 , (B.18) 

f d~v 8 
#11= j —P^A^A u P^A p A a = -A\ (B.19) 

r dv ~ 9 9 

#i 2 = / —P^A,A U V -AV-A = --A 2 A 2 + -A 4 . (B.20) 
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